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This paper derives Langevin equations for an optical cavity where the dominant nonlinearity 
arises from free-carrier dispersion. We define a generalized Wigner function, compute a Fokker- 
Planck equation that approximates the master equation, and convert this to a system of stochastic 
differential equations. These equations are similar to the Wigner equations for an optical Kerr 
cavity, but have additional noise terms due to the incoherent carrier excitation and decay processes. 

We use these equations to simulate a phase-sensitive amplifier and latch and compare the results to 
a Kerr model. 


Optical logic requires a platform that is fast, low-power 
and scalable to compete with electronics. In the past 
decade, nano-photonics has advanced to the point where 
optical cavities of size < (A/n)^ and Q factors > 10^ can 
be fabricated with standard techniques [1]. The hope is 
that these cavities can be used to amplify the optical non¬ 
linearity of materials or defects and perform all-optical 
logic for communications and computing at speeds and 
energy scales comparable to electronics. 

Free-carrier dispersion is a promising nonlinearity for 
low-power optical logic. The effect arises in all semi¬ 
conductors. In a semiconductor, there is a filled valence 
band and an empty conduction band, and when pho¬ 
tons are absorbed, they excite electrons from the valence 
band to the conduction band. Each absorption creates 
two free carriers - an electron and a hole - which evolve 
independently and decay on some timescale set by the 
material and its geometry. The carriers provide feed¬ 
back to the optical field by altering the absorption of the 
material (free-carrier absorption) or its refractive index 
(free-carrier dispersion). On timescales long compared 
to the free-carrier lifetime, it acts as an effective opti¬ 
cal nonlinearity and can be used to construct switches, 
amplifiers and other logic elements. 

Accurate, semiclassical models for free-carrier effects 
already exist, and these are valid when the carrier and 
photon number are very large [llEl- However, the real 
promise of free-carrier effects lies in their application to 
low-power photonic computing. In some materials, free- 
carrier effects are strong enough that switching can be 
achieved with as few as 100 photons per cavity. In this 
regime, quantum effects become important and place fun¬ 
damental limits on device performance. For example, 
quantum fluctuations in the photon number add noise to 
quantum amplifiers [4] and lead to spontaneous switch¬ 
ing in optical memories 016 ]. This motivates the need to 
develop a quantum model for the free-carrier nonlinearity 
that works at low photon numbers, similar to the mod¬ 
els that exist for cavity quantum electrodynamics (QED) 
and (Kerr) systems [71[8]. 

In this paper, we derive a quantum-mechanical model 
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for the free-carrier nonlinearity, following the standard 
open quantum systems formalism used for cavity QED, 
optical parametric oscillators (OPO’s) and systems 
[9]. However, simulating even a single cavity in this 
model is not practical, since the large number of available 
carrier modes makes the full Hilbert space exponentially 
large. Using a method based on the Wigner function 
UMa, we can reduce the master equation to a set of c- 
number Langevin equations that are simple to simulate. 
These equations bear resemblance to semiconductor laser 
rate equations and Bloch equations found in the litera¬ 
ture HMSj. An adiabatic elimination reduces the model 
further, giving a set of stochastic differential equations 
(SDE’s) for the field, electron number, and hole number 
in the cavity. The deterministic part of these equations 
matches the classical models found in the previous liter¬ 
ature, but the noise terms are new - and have a quantum 
origin. 

Section [T| introduces the quantum model for the free- 
carrier cavity. In Section [Hj we introduce the Wigner 
formalism and apply it to this model, deriving a set of 
stochastic differential equations (SDE’s) which we sim¬ 
plify by invoking a weak-doping, fast-dephasing limit. 
(Extensions of this model, incorporating two-photon ab¬ 
sorption and free-carrier absorption, are treated in Ap¬ 
pendix . 

The free carrier model we derive, summarized in Equa¬ 
tions ( [4Qp^ and ( 48p9 ), resembles the equations of mo¬ 
tion for the Kerr cavity derived in m, but there are 
extra noise terms. The steady-state behavior, and corre¬ 


spondence to the Kerr cavity, are treated in Section HI 


Next, the free-carrier SDE’s are applied to simulate two 
devices: a phase-sensitive amplifier in Section |IV| and an 
all-optical SR-latch in Section [V| Eor the amplifier, the 
free-carrier device does not show squeezing, whereas its 
Kerr analog does. Eor the latch, the spontaneous switch¬ 
ing rate is larger for the free-carrier device, and the dis¬ 
crepancy grows as the latch’s bistable states become more 
widely separated. 


I. QUANTUM MODEL 

Consider a single-mode optical cavity fabricated from 
an undoped semiconductor. The optical degree of free- 
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FIG. 1: (Color online) Electronic and optical modes of a semi¬ 
conductor cavity. 


dom can be represented as a harmonic oscillator, with the 
creation / annihilation operators a,a^. In band theory, 
the electronic degree of freedom is represented by many 
uncoupled fermion modes. For a two-band model, we 
have a single electron band and hole band. Each mode 
has its own (fermionic) creation / annihilation operators 
- Ck, for electrons, Ck, for holes, where k is the mode 
index. The operator algebra is: 

[a, a’l'] = 1, {ck, cj} = {ck, cj} = (5ki (1) 

The operator algebra is key to the Wigner analysis 
in Section |n[ In short, one can define a generalized 
Wigner function for the quantum system if one can find a 
“closed” set of operators {X^}, where the bath-averaged 
time derivatives ((iX^/dt)ad, {d{XiXj)/dt)Q,^^ defined in 
Section [nj are always functions of the {W}- However, 
Wigner functions for fermionic operators require the use 
of Grassmann variables m, for which the analogy to 
classical phase space is less intuitive. Thus, we identify 
fermion pairs and perform the following bosonization: 


CT-k 

— CkCk 

(2) 

C^+k 

II 

(3) 

77/k 

= 4 ck 

(4) 

rik 

= cj^Ck 

(5) 

Qk 

= UkUk = Cr+k<^—k — cj^cj^CkCk 

(6) 


This is similar to the operator algebra in an ensemble 
of two-level atoms, but there are some extra terms. Anal¬ 
ogous to the atom ensemble, the electronic polarization is 
given by cr_k. However, the free-carrier system contains 


two number operators rik, hk rather than one, as well as 
a pairing operator Qk- These arise because the electrons 
and holes in the free-carrier system have more freedom 
of movement: in an ensemble of atoms, each electron is 
confined to its parent atom and Uk = hk = Qk, while in a 
semiconductor these three quantities are no longer equal, 
since electrons and holes freely scatter between modes k. 
Operators ill) are bosonic because they are products of 
an even number of fermionic operators. Note that sev¬ 
eral bosonic operators, namely Ckcj^ and Ckc|., are not 
included in - this is because we are interested in 

systems that respect charge conservation, while Ckcj^ and 
dk^k violate it. 

The bosonized operators are closed under commuta¬ 
tion, with the following nonzero commutators: 


[cr±k,^l] = T<^kicr±k (7) 

[cr±k,hi] = T<^kicr±k (8) 

[c^±k, cr^l] = ±^kl(^k + f^k - 1) (9) 

[cr±k,Ql] = T^^klCr±k (10) 


Note that these are all commutators, rather than 
anti-commutators, because the operators have been 
bosonized. 


A. Hamiltonian 

The Hamiltonian consists an an optical part i7ph which 
resembles a harmonic oscillator, an electronic part Hq\ 
given by the electronic band structure, and an interaction 
part i^int due to the A • p light-matter interaction. It can 
be written as: 

H = AcO^a + I^Ak + ig^{a^(j_y. - aa+y,) 

V—^ K ^ ‘v*" ^ "V^" 

^ph ^el ^int 

( 11 ) 

where Ac = ^(co’c —cjph) is the cavity resonance detuning, 
Ak = E\^^c — E\^,v — ^ph is the detuning of the transition, 
and ^k is the atom-photon coupling. The coupling can 
be expressed in terms of material parameters as shown 
in Appendix [P] 

Hamiltonian ( pT| resembles the cavity QED Hamilto¬ 
nian. This is because both systems contain an optical 
term, and electronic term, and a light-matter interaction 
of the A-p form. Thus, it should not be surprising if free- 
carrier cavities exhibit many of the same phenomena ob¬ 
served in cavity QED, e.g. bistability, amplification, limit 
cycles [18]. 

B. External Interactions 

In the cavity, the optical field is relatively well isolated 
from its environment. The two primary interactions are 
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FIG. 2: (Color online) Major free-carrier effects in an optical 
cavity 


optical absorption, which gives rise to particle-hole pairs 
and is treated through Eq. 0, and coupling to the ex¬ 
ternal waveguide. Because these couplings are usually 
quite weak, the optical field tends to retain its coherence 
in spite of them. 

The same is not true for the carriers. Many forces act 
to dephase, thermalize, and scatter the free carriers on 
very quick timescales (typically around ^ 10 — 100 fs) 
Einaiio]. Even for very poor cavities with Q < 1000, 
this is much faster than the photon lifetime. The practi¬ 
cal upshot of this will be that, on optical timescales, the 
“coherent” part to the carrier fields cr±k can be adiabat- 
ically eliminated and only the “slowly-varying” carrier 
numbers Uk, ffk, Qk remain relevant to the system. 

Eor a bosonic, Markovian bath, external interactions 
can be treated by adding extra Lindblad terms to the 
Master equation [9]. The main external processes are 
given in Eigurej^ above. As Lindblad terms, they are: 


C. Single-Mode Theory 

To make the computation more tractable, we assume 
for now that all carrier modes k are identical. This is not 
always a realistic assumption, and in Appendix |C 1| we 
extend the result to non-identical modes. But assuming 
identical modes for now, Ak,5'k, and all the 7 k’s become 
independent of k. Suppose that there are N modes. One 
can now define mode-summed operators: 

a± = ^(T±k, etc. (16) 

k 

This reduces the dimensionality of the state space from 
bN -|- 2 to 7. In terms of these, the Hamiltonian and 
interaction terms are: 


H 


Acav 

Arc,k 

-Lnr,k 


L 


sCjk^l 


AcO^a + -Ae(n + n) + ig{a^a- 

- acr+)(17) 

a 

(18) 

\l Ire k 

(19) 

\/ Tnr Ck, \/^nr ^k 

(20) 

/ Tsc f / Tsc _f_ 

Y Y 

(21) 


Note that the sum operators used are bosonic, not 
fermionic. The single-mode theory would not work if 
we had started with the fermionic operators. 


II. WIGNER FUNCTION AND SDE’S 
A. Derivation from Quantum Model 


• Cavity Loss, mediated by 

Acav — \/^ ^ 


Under the quantum model described in Section [I| 
above, the state of the cavity is given by a density matrix 
p and evolves according to the master equation: 


• Recombination, mediated by 

Arc,k — \/7rc,k k (^3) 

• Nonradiative Decay, mediated by 

-Lnr,k — \/7nr,k ^k, ^ 7nr,k ^k (1"^) 

• Scattering / Dephasing, mediated by 

Ascjk^l — \/Tk^l Ck, \/7k^l cjck (1^) 

Some materials also have significant free-carrier ab¬ 
sorption and two-photon absorption. For simplicity, 
these are not treated presently, but are discussed in Ap¬ 
pendix |C 2[ 


^ = -i[H,p] + I ^ (2L,p4 - pL^p - pLllPj 

k 

( 22 ) 

Unfortunately, with an exponentially large Hilbert 
space, it is not practical to compute p or its evolution. 
To get around this problem, we express p in terms of a 
generalized Wigner function, compute the equations and 
show that they can be approximated by a Fokker-Planck 
equation, and solve the Fokker-Planck equation stochas¬ 
tically using SDE’s. 

This approach was derived by Carter [21] for optical 
fibers with a interaction; the same approach can 
be applied to optical cavities or cavity-based photonic 
circuits m- In both cases, there is an injective lin¬ 
ear mapping between optical states p and Wigner func¬ 
tions W{a,a*). Gronchi and Lugiato [11] extended this 
method to weakly-coupled many-atom cavity QED. In 
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this case, in addition to an optical degree of freedom, 
one must also keep track of N two-level atoms. The 
procedure is to define a closed Lie algebra {Xi} of both 
optical and atomic operators, and a semiclassical phase- 
space that with c-number variables {xi}. The generalized 
Wigner function is defined analogous to the optical func¬ 
tion - the Fourier transform of a characteristic function: 


W{x,t)= (23) 

In our case, the operator algebra consists of the optical 
and single-mode carrier operators (16) and is given as 
follows: 


X = [a, a’*', cr_, (7+, n, n, Q] (24) 

This is a 7-dimensional, operator-valued vector. The 
Wigner function thus lives is a 7-dimensional phase space, 
defined over the c-number variables: 


X = [a, a*, i;, i;*, m, m, q] (25) 

Under certain closedness conditions discussed in Ap¬ 
pendix the Wigner function satisfies a generalized 
Fokker-Planck equation. Truncating at second order, 
this reduces to a Fokker-Planck equation. The validity 
of this truncation generally requires that nonlinear cou¬ 
pling constants be weak m, and for two-level systems, 
that the number of atoms (carrier modes in this case) be 
large m- Both assumptions are true here. As a solution 
to a Fokker-Planck equation, the Wigner function can be 
sampled stochastically by solving the following SDK’s: 


Computing the cumulant matrices and from 
the quantum model is straightforward but very 

tedious, so we used Mathematica to derive the result. The 
details are given in Appendix [B] the SDE’s are: 


da 

dv 


dv^ 


dm 

dffi 

dq 


- iAc'j a + gv dt + d^i (30) 

(^-^ + iAc'J a* + gv* dt + d ^2 (31) 

—a{N — m — ffi)g 

(^-1^ - iA.'j v\dt + d^3 (32) 

—a'^{N — m — ffi)g 

iAe) V* dt + d^i (33) 


+ 


+ 


(-^+iAej V*'^dt + dU 

[-g{av"' -h va^) - ^nrm - ^rcOl dt + (34) 

[-giav'^ -h va'^) - ^nrffi - -frcOl dt + d^6(35) 
—1sc{q — mfh/N)dt + d^j (36) 


where jtot = 7sc + 7rc + 7nr + 7nr and the noise processes 
d^i have the covariance matrix: 


d^d^^ = dt 


(37) 


where Ccav + + Crc^ + Cn?^ + is the sum of 


the terms in Eqs. (B^B6). Note how Equations (30 


36) resemble both the Maxwell-Bloch equations and the 


Langevin equations for many-atom cavity QED derived 
by Gronchi and Lugiato m- However, because of scat¬ 
tering between carrier modes, we need to keep track of 
m, fh and q separately. 


dxjYi — gm dt -\- ^ ^ RfYin dw^iq 


(26) 


B. Approximations 


with dwn a Wiener process and 


1. Fast Dephasing, Nondegenerate Exeitation 



Three approximations make these equations more 
tractable: fast-dephasing, nondegenerate excitation and 
the single-carrier approximation. Fast dephasing assumes 
that the scattering rate jsc and detuning Ae are faster 
than any other timescale in the system, thus 

'Inrt 'Irc’) '1 fc 7sc: (38) 


where the time-derivatives are computed using the ad¬ 
joint to (22) 


= -i[^, 2LlALk - LlLkA - 

k 

(29) 

and (•••)?? is defined so that normally ordered prod¬ 
ucts return simple c-number polynomials, e.g. (a)p = a, 
^{o)a + ao))p = (a*(a, {o)a)p = a^'a — etc. 


dA 

dt 


This is related to the relaxation-time approximation 
that holds for most semiconductors [2]. In useful, opti¬ 
mized free-carrier devices, all of the carrier timescales 
“ Inr^lrc^lfc “ are of order the photon lifetime. To 
achieve strong carrier effects, we generally have cavities 
with Q > 1000, giving a photon lifetime of Tp}i > ps. 
Ultrafast studies show that inter-mode scattering takes 
place on times of order 10-100 fs [191 ED], giving scat¬ 
tering rates at least 10-100 times faster than any other 
timescale in the system. 
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Fast dephasing leads to an adiabatic elimination of the 
dipole terms and the pair density q. These vari¬ 

ables will be replaced by their steady-state values, and 
a new set of SDE’s are obtained for the reduced basis 
(a, a*, m, fh). 

Nondegenerate exeitation assumes that the number of 
carriers is much less than the number of carrier modes 
- in other words, the valence and conduction bands are 
far from being degenerately filled with electrons or holes. 
This approximation is invoked by setting 


m,m,q N (39) 

This simplifies the equations of motion by discarding 
effects like absorption saturation that are negligible for 
low-power, high-Q optical logic devices. The resulting 
equations of motion are: 


da 

dm 

dffi 


2 


i(Ac + 5 cm + 5 cm) 


adt + d^a 


[rja'^a — 7nr^ ~ dt + d^rn 

[rja'^a — — 7 rc^^] dt + 


(40) 

(41) 

(42) 


with noise terms 


d^a = -V^d^in - (43) 

TTi 2-^/7 Fte [or H- '\Jdw^yi 

dWrc (44) 

d^fh = 2^Re[a*(i/3^] + y/jnr^^dWfn 

T^y^c^^ dWrc (45) 


where the d/d’s are complex vacuum Wiener processes, 
e.g. d/3*d/d^ = and the dre’s are real Wiener pro¬ 
cesses, e.g. dw‘^ = dt. 

In the equations above, we rescaled y^c and defined a 
bandfilling carrier-dependent detuning 5c and linear ab¬ 
sorption y: 



T] = 2ATm[dc] (47) 


Since Ae = cjg — cj, this function has one pole (for uj) 
in the lower half-plane, uj = uje — As a result, its 

real and imaginary parts satisfy the Kramers-Kronig re¬ 
lations. The carrier-dependent dispersion and absorption 
are given by the real and imaginary parts of dc, respec¬ 
tively. Around l/5c carriers are needed to shift the cavity 
resonance by one linewidth; since this quantity is much 
smaller than N under the nondegenerate approximation, 
it follows that Im[dc] ^ Re[dc]. For a pure bandfilling 
effect, we can generally neglect the imaginary part. 

The d/din in ( [4Q][4^ is the vacuum noise of the input 
field and d/3^ is the noise due to linear absorption; each 


behaves as a vacuum Wiener process d/d* d/3 = |dt m- 
The dwm^ dwfh and dwq are real-valued noises due to 
carrier loss and recombination, and go as dw^ = dt. 

The noise term for a is fairly standard for open quan¬ 
tum systems: a sum of two vacuum noises. The noise 
terms for m and fh have Poisson statistics: for each pro¬ 
cess with rate Rdt^ there is a corresponding noise term 
\fRdw. Since carrier generation involves photon absorp¬ 
tion, one should not be surprised by the Poisson noise on 
this signal. Likewise, since the carrier number is quan¬ 
tized and carrier decay is a random process, there should 
also be Poisson noise on the decay terms. 


2. Single-Carrier Approximation 


In many cases, the equations (40 [^) can be reduced 
further by positing that m = fh and introducing an ef- 
feetive earner number Nc equal to this quantity. For 
example, it will hold if only one of the carrier species is 
relevant (for instance in silicon, where Sc » Sc due to 
the plasma effect 0 ), if the recombination process y^c 
is dominant, or if the number of recombination sites is 
limited (much smaller than the number of carriers) so 
that electrons and holes tend to decay together. If any 
of these cases hold true, equations ( 4Qp2 ) become: 


da 

dNc 


K + T] 
2 


^(Ac + 5cNc) 


a dt d^a 


(48) 


[7]a*a- JnrNc - -frcN^] dt + d^N (49) 


with noise terms 


d^a 

= -VKdpic,- .Jrjdp,^ 

(50) 

d^N 

= 2y^Re[Qf d/Sp] T '\J'^jij-NcdWjic 



+ \/'1rcN'^dWrc 

(51) 


In the sections below, we consider the behavior of op¬ 
tical logic devices assuming the single-carrier approxima¬ 
tion. Note that this is not strictly necessary - it is equally 
feasible to use Eqs. ( 4Q]|^ ) to obtain similar qualitative 
results. 


III. STEADY-STATE BEHAVIOR 


Consider now the case where ( 48][^ ) hold. Suppose 
also that direct recombination is negligible {jrc = 0)- 
This is the limit to which III-V photonic crystals op¬ 
erated near the band edge, which have the best perfor¬ 
mance to date m, belong (the parameters in Table are 
less than a factor of 5 from the state of the art). In this 
limit, a free-carrier cavity modeled by equations ( 48p^ 
has two timescales - an optical lifetime Tph = 1/ (^ + y 
and a free-carrier lifetime Tc = 1/ynr- 
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Tunable 
Phase Shifter 


Cavity with x^^^-like 
nonlinearity 


\ 





a 



1 real degree of 
freedom: carriers 


1 complex degree of 
freedom: optical field a 


FIG. 3: (Color online) Adiabatic elimination of a free-carrier 
device into a tunable phase shifter (left) and a Kerr-like non¬ 
linear cavity (right). 


First, a steady-state limit is discussed. This is the 
case when both the carrier and photon lifetimes are much 
shorter than the relevant timescales. Questions of ther¬ 
mal stability, for instance, can be treated in the steady- 
state limit. Next the limiting case of Tph ^ Tc, where the 
carrier population varies much faster than the photon 
population, is treated and we show that the free-carrier 
model reduces to a Kerr model with extra noise terms. 

In this section, we work in normalized units by setting 
k = K T] ^ 1. Rates, time constants and powers are 
scaled by appropriate powers of k. This allows our results 
to generalize to a wide range of systems spanning orders 
of magnitude in speed and size. 


A. Steady-State Limit 

In the steady-state case, we set all noise terms to zero 
and solve for TVc = d = 0. Solving for d = 0, the steady- 
state internal field a can be related to Nc and the input 
field Pin as follows: 


-y^lSin 

+ i{A(. + ScNc) 


(52) 


This is the familiar formula for the field in a resonant 
cavity, where the detuning Ac + ScNc, depends on the 
free-carrier number. Solving the Nc = 0 equation gives 
This can be rearranged into a polyno¬ 
mial equation for a*a: 


= {a*a) {k/2f + (Ae + {vSc/lnr){a*a)y 


( 53 ) 


Name 

Description 

Value (this paper) 

k 

Total linear loss k = K-\-r] = uj/Q 

0.5 ps-^ 

K, 

Output coupling 

0.8k 

r] 

Linear Absorption 

0.2k 

P 

2-photon Absorption 

0 

X 

Kerr 

0 

Sc 

Free-Carrier Dispersion 

0.013k 

Xrc 

Recombination 

0 

Xnr 

Free-Carrier Decay 

1.0k 


TABLE I: System parameters used in the free-carrier simula¬ 
tions in this section. Values given are similar to state of the 
art III-V photonic crystals. 


When the external power P = Pi^Pin is set, this is 
a cubic equation for the internal photon number = 
a*a. It is the same optical bistability cubic as the Kerr 
cavity mill], with the effective Kerr nonlinearity: 


Xeff — 


Inr 


(54) 


One can solve the cubic (53) to obtain a*a; it is not 


always uniquely defined. Just like Kerr cavities and atom 
cavities, free-carrier cavities exhibit hysteresis and bista¬ 
bility, with both “low” and “high” intensity states being 
allowed for the same input power. Figure shows the 
stable lower- and upper states, and an unstable middle- 
state, for varying values of Ac- 

The intuition behind this bistability is that, when the 
cavity is off resonance and a sufficiently large number of 
carriers are injected, it will shift back on resonance. If 
there is a strong enough input, then a large power builds 
up inside the cavity and this large carrier population can 
be maintained through absorption, giving rise to the high 
state. On the other hand, if there are no carriers to begin 
with, the cavity stays off resonance and there is never 
enough power in the cavity to raise the carrier number - 
hence the low state. Analytically, one can show that the 
bifurcation sets in when: 


Ac <-^1 (« + ??) (55) 

Much of our intuition behind free-carrier nonlineari¬ 
ties comes from this steady-state picture. It does not in¬ 
clude any quantum effects or even any dynamics, but the 
shapes of the curves in Figure [^suggest that the device 
could be used as an amplifier or a switch. We will show 
in a later paper that free-carrier cavities can do much 
more than this, but that will build on the fundamentals 
discussed here. 

The steady-state picture has been amply discussed in 
the literature miiHiiii], so it is not worth describing in 
more detail here. Rather, we now proceed to look at the 
quantum noise and dynamics of these systems. 
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FIG. 4: (Color online) Steady-state solutions to Nc for optical 
free-carrier cavity at different detunings. 


B. Effective Model 


Next, we go from the steady-state picture to the limit 
of short carrier lifetime. In this opposite limit, Tph ^ 
Tc- Typical devices do not realize this limit, but it is 
useful because it enables an apples-to-apples comparison 
between the free-carrier and Kerr effects. 

To adiabatically eliminate the carrier number, one re¬ 
places Nc with its steady-state value: 


'Inr 


(56) 


This gives the following SDE for the relevant dynamical 
variable, a\ 


da = 


i f Ac ff a dt-y/^dPin+d^'^ 

^ \ '^nr J 

(57) 

where the is a new noise term that depends both on 
the d^cK and As before, the analogy to the Kerr 

model is clear: Equation (57) is very close to the Wigner 
equations for the Kerr cavity m, but the noise term is 
different. The effect of this noise term will be discussed 
in the following sections, where the performance of Kerr- 
and free-carrier based amplifiers and switches is analyzed. 



Outl (Signal) Ini (Signal) 



FIG. 5: (Color online) Left: Simple single-cavity amplifier, a 
cavity sandwiched between phase shifters and displacements 
G — L{^') < < (Cav) < < L{I3). Right: Symmetric two- 

cavity amplifier. 


IV. AMPLIFIER 


Figureshows that, for certain detunings, the state of 
the cavity changes very rapidly with a change in input 
power. One can imagine using such a device to amplify 
differential signals: if the input signal is perturbed, that 
perturbation will be multiplied by some gain factor in 
the output. 

The real picture is actually a bit more complicated, 
since the input field has two quadratures. In the Kerr 
cavity, one of the quadratures is amplified while the other 
is de-amplified [24]. This gives rise to phase-sensitive 
amplifieation which, since there is no additional noise in 
the Kerr system, also squeezes the quantum noise of one 
quadrature below the vacuum level. 

Key to an optical amplifier are its gain its noise 

spectrum S{uj), and the scale on which nonlinear effects 
take over. The gain and noise can be predicted by lin¬ 
earizing equations of motion ( 48p9 ) around the steady- 
state value. This takes the general form: 


dx = Ax dt ^ B d/^in ^ F dw (58) 

dpout = Cx dt F D dPin (59) 

where x and P are doubled-up state vectors, which include 
the complex field operators and their conjugates [25], as 
well as the (real) carrier number: x = ((5(a, 6a*, 6Nc), 
dp = {dp, dp*) (removing any constant coherent input), 
and a, Nc are the steady-state values. 

Linearization is key because many general results of 
stochastic systems theory only apply to linear or approx¬ 
imately linear systems [26]. For example, in a linearized 
system, the output squeezing spectrum can be computed 
exactly for Gaussian inputs [22] ED. Many results in 
quantum feedback control theory are also restricted to 
linear systems [28] [29] . 

With a linearized model in hand, it is a simple matter 
to compute the internal state covariance a, the transfer 
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FIG. 6: (Color online) Plot of the maximum and minimum 
gain G+(ct;), G-(uj) for free-carrier cavity (left) and Kerr cav¬ 
ity (right). Parameters are from Tablewith Ac = —0.7. 


and noise matrix T{uj)^ N{uj)^ and the frequency-domain 
input-output relation ESHSOj: 


free-carrier cavity becomes very efficient at amplifying 
off-resonance, whereas the Kerr cavity hardly amplifies 
at all. 

Gain is maximized when the system is very close to 
instability - that is, when at least one of the eigenvalues 
of A is very close to the imaginary axis. From (61), an 
eigenvalue decomposition of A gives the transfer function 
the following form: 


t = d + Y. 

i 


—iuj — Xi 


(62) 


where Vi^Ui are related to C and the eigenvectors and 
are the eigenvalues of A. Since B,C ^ O(y^), the 
numerator term is proportional to k. Near the resonance, 
the sum is dominated by the eigenvalue closest to zero, 
Amax- The maximum gain should intuitively take the 
form: 


G- 


0{i^) 

I -^max I 


(63) 


This is a Lorentzian with a peak at Im(Amax) and 
bandwidth of Acj = —Re(Aniax)- The peak gain is 
thus Gmax = —0(/i:)/(Re(Amax))- This gives us a gain- 
bandwidth relation: 


G^axAcj = 0(/i:) (64) 


Aa + + FF^ = 0 


Pout,UJ 


D^C— 


1 


iuj — A 


tB 


An,a; + G- 


1 




-iuj — A 

N(uj) 


(60) 


(61) 


Unfortunately, because the doubled-up matrices here 
are 3-by-3 rather than 2-by-2, the analytic results are 
rather cumbersome and therefore not reproduced here. 
Instead, in this section I compute these quantities numer¬ 
ically and compare the results to the Kerr system. The 
results here are compared against a Kerr cavity with the 
same effective nonlinearity, x = 


Gain 


The gain is computed from the singular values of the 
doubled-up transfer function T{uj). If both singular val¬ 
ues are the same, the device is a phase-insensitive am¬ 
plifier. Both the Kerr and free-carrier cavities, however, 
only amplify one quadrature. As Figure [TV A shows, they 
de-amplify the other quadrature as well. 

At and below the ideal input ~ 9, the Kerr and 
free-carrier cavities seem to amplify in the same way. For 
over-driven cavities, the behavior is very different. The 


The greater the amplifier gain, the slower it responds 
and the narrower its bandwidth. 


B. Internal State 


The internal state is computed using the Lyapunov 
equation (60). This time, the Wigner equations contain 


additional noise terms, which make the state noisier than 
the state of an equivalent Kerr cavity. This is plotted in 
the Figures [7][^ The state remains roughly Gaussian, 
but the size of the Gaussian is larger than in the Kerr 
case, especially above the inflection point. 

Unlike in the Kerr case, the mode in the free-carrier 
cavity is never squeezed. As seen in Figure the eigen¬ 
values cr+, cr_ of the covariance matrix a are always > 
ensuring that the state is always “classical” in the sense 
that it has a valid P representation. Given that the car¬ 
rier excitation and decay process is highly incoherent, it 
should not be too surprising that the cavity always re¬ 
mains in a classical state. But it is a clear departure from 
the Kerr model, and this classicality could conceivably be 
used to distinguish between the two in an experiment. 

Also note that the noise grows linearly with the in¬ 
put field at high powers. This happens because the free- 
carrier number is constantly fluctuating, being driven by 
excitation and decay events that mimic a Poisson pro¬ 
cess. At high carrier numbers, this means that the cavity 
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FIG. 7: (Color online) Internal state of free-carrier cavity, 
simulated eigenvalues of a ((7+ and cr_, the larger and smaller 
eigenvalue, respectively) compared to analytic result (solid 
lines). The dashed line is the analytic result for an equivalent 
Kerr cavity. Ac = —0.7 



ijj (normalized) 


FIG. 9: (Color online) Noise spectrum modes S- for 
the a free-carrier cavity with Ac = —0.7 at various pump 
powers. Green and blue lines are numerical simulations; red 
solid line is the prediction from the linearized ABCD model. 
The dashed lines are the prediction from the Kerr model. 



An =5.00 


An =9.00 An = 13.00 An = 17.00 


- 2-10 1 2 - 2-10 1 2 - 2-10 1 2 
Re[a-ao] 


FIG. 8: (Color online) Simulated Wigner functions for free- 
carrier (top) and Kerr (bottom) cavities with the same effec¬ 
tive - Analytic approximation for linearized model given 
in white contours. Ac = —0.7 


detuning and consequently the cavity field become very 
noisy. This does not happen in the Kerr cavity, where the 
nonlinearity is mediated by virtual transitions which do 
not add any noise to the system. It is a peculiar conse¬ 
quence of the incoherence of the free-carrier mechanism. 


Carriers Kerr 



C. Output Noise Spectrum 


FIG. 10: (Color online) Noise spectrum *S'+, S- as a function 
of input /3in and frequency u for Kerr and free-carrier models. 
Ac = -0.7 


Given a linearized input-output model, we can com¬ 
pute the squeezing spectrum (noise spectrum) for the 
cavity output field [laiii]. The squeezing spectrum for 
quadrature 0 is defined as the power spectral density of 
a homodyne measurement of /5out(^)- That is, for the 
following homo dyne signal. 


(65) 


the squeezing spectrum is: 

S0{lo) = ./2P^), ( 66 ) 

27rd[uj — uj') 

Se{uj) is normalized so that the coherent state has 
Se{oj) = 1. For general states, Se{oj) depends on 0. The 
maximum and minimum of So{uj), with respect to 0, are 


jeit) = e 
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FIG. 11: (Color online) Circuit diagram for a photonic SR- 
Latch 


denoted S-^{uj) and S-{uj)^ respectively. 

The squeezing spectrum of the Kerr cavity can be com¬ 
puted analytically m- By contrast, since the free-carrier 
squeezing spectrum involves the inverse of a 3 x 3 ma¬ 
trix, it is unlikely that a simple expression can be found. 
However, it is not difficult to compute numerically. 

In Figure the noise spectrum is obtained in two 
separate ways: first, simulating the full system in the 
time domain and taking the Fourier transform of the ho- 
modyned output (blue, green curves); and second, from 
the analytic predictions of the linearized ABCD model. 
These agree everywhere except for very large pump pow¬ 
ers, where the system approaches a bifurcation. 

Figure displays the noise spectrum for the whole 
range 0 < /3i^ < 20, for both Kerr and free-carrier de¬ 
vices. Two things are obvious. First, the noise curve (at 
least for the 5'+ component) matches the general form of 
the gain curve in Figure [TV A[ This is of course necessary 
because there must be noise wherever there is gain. The 
free-carrier cavity, unlike the Kerr cavity, amplifies not 
only at cj = 0 near the point of maximum gain, but also 
for cj 7 ^ 0 for above that point. 

Unlike the Kerr cavity, the free-carrier cavity does not 
squeeze the output field. Regardless of the parameters, 
regardless of the pump power, both 5'+ and S- are al¬ 
ways above the vacuum level, indicating that this is a 
classical field with no squeezing. The Kerr cavity, on 
the other hand, squeezes light over a broad range of the 
spectrum. This is in agreement with the results of the 
previous section, which showed that the internal field of 
the free-carrier cavity was classical. If the input and in¬ 
tracavity field are in a classical state, so is the output. 


V. LATCH 


It is also possible to construct a switching device us¬ 
ing only amplifiers, provided the amplification is large 
enough [32]. The circuit in Figure 11 uses two identical 


amplifiers in a feedback loop. Suppose that each ampli¬ 
fier has a gain G. Consider the fate of a perturbation in 
the top amplifier. An input SP is amplified to G Sp. This 
amplifier has a fan-out of 2, so {G/V2)6f3 passes to the 
right and exits the system, while {G/\/2)6P passes to the 



FIG. 12: (Color online) Top: time series of the latch internal 
state. Bottom: input field. 


lower amplifier. 

In the lower amplifier, it grows to (G^/2)^/3, is fed back 
into the original amplifier. After passing through this 
loop, the signal strength has grown to {G‘^/2)6/3. This 
leads to a latching instability if the gain is sufficiently 
large: 


G>V2 (67) 

Symmetry gives the latch some very desirable proper¬ 
ties. Unlike the single-cavity switch, the two states here 
are symmetric. Thus, there is less worry about finding 
the right bias field to “balance” the low and high state, 
and transitions between the states look the same. But 
this comes at the cost of the added complexity of two 
cavities, plus the extra connections. 

Figure [T^ shows a latch simulation for the same param¬ 
eters used in the previous section. Here, the detuning is 
set to Ac = —0.7, large enough to realize a large gain, 
but not large enough make an individual cavity bistable. 
The symmetry between the two states is very clear. 

Externally driven switching in the latch is good, be¬ 
cause it allows the user to set the state of the latch, which 
becomes a memory element. But thanks to quantum 
noise, Kerr and free-carrier devices also undergo sponta¬ 
neous switching. This is generally bad, because it limits 
the lifetime of a carrier-based memory. 

In the Kerr case, spontaneous switching is driven by 
vacuum fluctuations m- In the free-carrier case, vac¬ 
uum fluctuations combine with stochastic carrier excita¬ 
tion and decay to drive the switching process. Because 
there are more fluctuations, we naturally expect the free- 
carrier cavity to spontaneously switch at a higher rate 
than the Kerr cavity. 

When the switching rate is low, the switching process 
is well described by a two-state Markov chain. In a two- 
state Markov chain, there are two states a and b, with 
jump probabilities 


P{a b) = '^adt, P{b a) = ^^dt (68) 
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FIG. 13: (Color online) Asymmetric part of the latch state 
Re[ai — a 2 ] for free-carrier based latch (left) and Kerr- 
based latch of the same ■ Cavity detuning set to Ac = 
0.50,0.52,0.54. 



(normalized) 


FIG. 14: (Color online) Spontaneous switching lifetimes Tsw 
for the Kerr and free-carrier latch as a function of detuning. 


In the latch, the states are symmetric, so Xa = lb = 
'Xsw The probability of being in a given state evolves as: 

dPa_ p p 

^snjPa “t“ ^swPb 

^ = IswPa - IswPb (69) 

Solving this linear system, one finds that the system 
reverts to its equilibrium distribution with a character¬ 
istic time Tsw = 1/(27sw). This time can be measured 
from simulations of the latch by looking at the autocorre¬ 
lation function i?(r), which decays exponentially for the 
Markov process: 


to —0.54, about where the latching transition occurs. 
Larger negative detunings correspond to higher gain (see 
Fig.0, and likewise stronger latching. However, for a 
fixed detuning, the free-carrier cavity has a shorter spon¬ 
taneous switching lifetime. 

This is also seen in Figure which plots Tqw foi* the 
free-carrier and equivalent Kerr latches. Because the free- 
carrier cavity has more quantum noise than the Kerr cav¬ 
ity, its spontaneous switching rate is higher. The effect 
becomes noticeable once the latching transition sets in, 
and grows as the latching grows stronger. 


VI. CONCLUSION 


In this paper, we introduced a method to simulate opti¬ 
cal cavities where free-carrier dispersion is the dominant 
nonlinearity. This method is based on deriving an ap¬ 
proximate Fokker-Planck equation for the Wigner func¬ 
tion, the approximation being valid in the weak-coupling 
limit where the detuning per carrier is much smaller than 
the cavity linewidth and the mean photon number is 
large. Importantly, this allows us to keep track of the 
dominant quantum effects (vacuum noise in the optical 
field, Poisson noise in the carrier excitation and decay) 
without running a full quantum simulation. This method 
was then applied to simulate an optical amplifier and an 
SR-latch. 

Because the semiclassical properties of these devices 
are well known, in this paper we focused on the quan¬ 
tum noise in the free-carrier amplifier and latch. Since 
the free-carrier dispersion creates an effective non¬ 
linearity, in both cases the free-carrier simulations were 
compared to simulations of an analogous Kerr-based de¬ 
vice, for which the quantum model is well known. In the 
semiclassical, steady-state limit, the Kerr and free-carrier 
devices behave the same. 

For the phase-sensitive amplifier, we find that the free- 
carrier and Kerr devices have the same gain near the 
amplification maximum. However, the noise in the am¬ 
plification direction is larger for most frequencies. In 
the other quadrature, the Kerr cavity squeezes the field, 
while the free-carrier cavity does not produce any squeez¬ 
ing. In the SR-latch, we notice a difference in the the 
spontaneous switching rate: the free-carrier latch has 
switching rates of > 5 larger than the Kerr latch, a dis¬ 
crepancy that grows as the latching becomes stronger. 
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Appendix A: Closedness of Operator Algebra 
1. Closedness and the Wigner Function 

The single-mode model of Sec. |I C| reduces the num¬ 
ber of phase-space dimensions from bN -|- 2 to 7, but the 
operator algebra X = [a, a’*', cr+, cr_, n, n, Q] is a very re¬ 
stricted basis set. Many degrees of freedom cannot be 
expressed in terms of the Xi. However, if certain closed¬ 
ness conditions are satisfied, operators in the algebra stay 
in the algebra under time evolution. Since the Wigner 
function is tied to expectations of operator products, this 
allows us to set up a PDE for the Wigner function. In 
essence, the degrees of freedom contained in the single¬ 
mode Wigner function exactly “decouple” from the other 
degrees of freedom in the system, and the single-mode 
model is valid. 

X is closed under commutation and thus forms a valid 
basis for an algebra S - a vector space spanned by the Xi 
and their products, e.g. XiXj, XiXjXj^, etc. We say that 
B is closed under time evolution if the time derivative of 
every element of S is in S: 


BeB 


dB 

dt 


e B 

ad 


(Al) 


If ( |A 1 | ) holds, then for every c-number product 
Xm ... there exists a polynomial such that 


(A 2 ) 

and from the correspondence between Wigner moments 
and operator products, we obtain an equation of motion 
for the Wigner function’s moments: 


d 

dt 


{Xm • • • ^p)w — {^m...p{^))w 


(A3) 


It is a well-known result in stochastic calculus that we 


for W(x, t), w. 
lants. Eqs. (26 


can recast ( [A3 ) as a generalized Eokker-Planck equation 
lere the moments are replaced by cumu- 


to second order. 


28) arise when this equation is truncated 


Eor example, Ck G Ak, Q E B^^\ G B^^\ Below, we 
prove several lemmas about the ordering of bosonic and 
fermionic operators. 

Lemma 1. If / G Ak and b G B^^\ then fb = bf-\- f' 
and f^b = bf^ + (r)^ where /', f" G Ak- 

Proof. The case for / is proved by a search of all rele¬ 
vant cases. B and Ak have 7 and 10 basis vectors, respec¬ 
tively, so this is 70 commutators to check (most are zero). 
Given this, the case holds because [f^,b]^ = — 

Lemma 2. If / G Ak and b G B^'^\ then fb = bf P 
EiKfl and pb = bp + EibUpy, where G Jk 

and 6',6" G 

Proof. Induction on n. The n = 1 case is proved 
in Lemma I. Assuming it holds for n — I, write b = 
J2jbi,jbn-i,j with bij G B^'^\bn-ij G Using 

both Lemma I and the n — 1 case, we move the fermionic 
operator from the left to the right side of the expression 
(summation signs omitted for claity): 


fb = fbijbn-l,j = {bijf + fpbn-l,j 

= h,j{bn-l,jf + bP 2 ,jkfjk) + {bn-l,jfj” + K- 2 ,jkfjk) 

(A7) 

With appropriate index renaming, this takes the desired 
form. The p case is analogous. 

Lemma 3. If /k,/k G Ak, then /j[/^ G Sk and 

Proof. Done by a search of all relevant cases - 100 in 
all since Ak has 10 basis vectors. 

Theorem 1. The operator algebra B is closed under 


(29) if the Hamiltonian is in B and the Lindblad terms 


take the following form: L ^ bf\^^ ... /k^, with b E B and 
either /k, G or G Jici, or i ~ 65 k, iov b G B,g e 
Sk- There must be one L for each multi-index k^. 

Proof. To prove closedness, we must show that ( j^ 
is in B for all A E B. This is the sum of a Hamiltonian 
and Lindblad terms. The Hamiltonian term —i[A^H] is 
obvious since both A and H are in the algebra of A, 
which is closed under commutation. 

The Lindblad term is \{2L'^AL - L'^LA - AL'^L). We 
first use Lemma 2 to move the indexed parts /k^ ... /k^ 
to the same side of the expression; for instance, for L^ AL^ 
we find 


2. Proof of Closedness 

We will prove closedness for a relatively broad class 
of Hamiltonians and Lindblad terms. To start, define a 
boson space Sk, a restricted fermion space Ak, and sum- 
operator spaces B^'^^ (note no index k): 


^ ■ ■ -in 1 • • 'jn ^ ^j 1 f 


(AS) 


To each term in this sum, we apply Lemma 3 to combine 
the fermionic operators into bosonic operators. 


Bk = span(cr+k,o-_k,nk,hk,(5k) (A4) 

J^k = {a;iCk + a;24 I xi,X2 e Bk} (A5) 

= span(Aii... < n), B = B°° (A 6 ) 






(/l,kJ^--(/4.kJ^&i/i.,k..../i„.k„ (A9) 
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Summation over k is critical here; without it bi would not 
be a bosonic sum-operator in Thus, the algebra B 

is not closed for a Lindblad term with just a single k - we 
must sum over all the k’s in order to recover closedness. 

Now we use Lemma 3 to move the bosonic operator b 
to the left, recombine the operators with index k 2 , and 
repeat until all fermionic operators have been combined. 
This gets rid of all the indices k^, resulting in an operator 
that lives in B. The terms LA and AL^L are done 


the same way. It follows that the Lindblad term in (29) 
lives in B. As before, the action of a single Lindblad 
term breaks closedness, but when we sum over k, it is 
recovered. 

The result for L = bg\^ can be shown without Lemmas 
1-3. We just use the commutation relations of the X\^ to 
move the all the indexed terms to the same side, where 
they can be combined and summed into a term in B. 

This theorem encompasses all the quantum models 
studied in this paper. A few examples of things it 
does not apply to would be index-dependent effects, say 
H or certain effects that violate charge conser¬ 

vation, such as L ^ Ckcj^. 


J 


Appendix B: Full Wigner SDE’s 

The cumulants and are computed from Eqs. 
( 27p8 ) using Mathematica. The terms are separated by 
physical origin in the sections below. 


1. Uncoupled Cavity, Carrier Terms 

In this case, H = AcO^a + Ae(n + n)/2 and L + ^/Ka. 
It is easy to show that: 


= 




'0 

K 

2 

0 

0 

0 

0 

o' 

(iAe - f )a* 


K 

2 

0 

0 

0 

0 

0 

0 

—iAeV 


0 

0 

0 

0 

0 

0 

0 

iAeV 

, = 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 


(Bl) 


2. Photon-Carrier Interaction 


Here, H = ig{a^(j- — a( 7 +), and there are no environment couplings. There is no noise term here. 


= 


gv 

gv^ 

—ga{N — m — ffi) 
—ga*{N — m — fh) , 
—g{av'^ -h ua*) 
—g{av* -h va*) 
—g{av'^ -h ua*) 


= 0 


(B2) 


3. Pree-Carrier Dispersion / Absorption 

Here, Lk = ^Jljca ^k <^5 -\Jlfca ^k^^ and H = a'*'a((5jcd^ + ^fcd^)- Define Sfc = Sfcd — Then the cumulants 

become: 


C(i) 


—i{Sfcm + Sfcfh)a 

-i{Sfc + Sfc)(a*a - \)v 

0 

0 

0 


C(2) 


0 (7/ca«^ + 7/ca«^)/2 0 (7/ ca T 0 0 0 

* 0 (jfca + 7fca)va* 0 0 0 0 

0 >K 0 (Vca + 7 /ca)(a a-^2)9 Q 0 

X (A -f 2g — m — m) 


* 0 

0 0 

0 0 

0 0 


* 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 


(B3) 


Equation (B3) has a nontrivial noise matrix. However, 
dephasing limit usually taken. 


this can be greatly simplified in the non-degenerate, fast- 
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4. Recombination 

This is mediated by the term L = cr_. Recombination only takes place when an electron and hole occupy the 
same state k, so the rate goes as the pair density g, not as the carrier density m + m. 


= 


Trc 


0 


'0 

0 

0 

0 

0 

0 

0 ' 

0 


0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

+ 2g — m — m) 

b 




, = 

0 

0 

l{N^2q- 

m — ffi) 0 




-Q 


0 

0 

7 

-7;* 

2 ^ 

Q 

Q 

Q 

-Q 


0 

0 

7 

2 ^ 

Q 

Q 

Q 

. 


0 

0 

7 

^7;* 

2 ^ 

Q 

Q 

Q _ 


(B4) 


5. Nonradiative Decay / Excitation 

Nonradiative decay is mediated through a term of the form L = yjlnr c/e • Strictly speaking, one must 

include write L = etc. where ri is the electronic mode into which the carrier decays, to make the L operator 

bosonic. However, if there are many more recombination sites than carriers, this mode’s dynamics are not relevant 
and the fermionic L gives the right result. 






'0 

0 

0 

0 

0 

0 
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0 
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0 

0 

0 

2 (7nr T 

, C(2) = 

0 

0 

0 

^InriN - ffl) 

+ linriN - m) 


2'l'rlr^ 

2 {inr T 

2 {.Unr T 

-Inrm 

0 

0 

2Tnr(-^ 

+ \lnr{N - m) 

0 

^InrV* 


2 i.'ynr T 'Inr)'^ 



0 

0 
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0 

'InrQ. 
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0 
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\lnrV 

\'inrV* 

0 

Inrrn 

7nr^ 
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2 

'InrQ. 

'JnrQ 

(7nr T 'Jnr^Q _ 


(B5) 


Scattering 


The scattering terms are L = c\^c\. Defining an average scattering rate by + 7sc) ^ 

75 c, we have: 


(7(1) 


7sc 


0 


'0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 



0 

0 

J_7;2 

|7Vei 

0 

0 


— -7;* 

2 ^ 

0 

0 

iTVei 

J_7;2 

0 

0 

lv^e2 
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0 

0 

0 
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0 
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Q + 


(B6) 


where ei = 1 — (m. + m)/N + 2 mm /62 = 1 + {2q — m — m — 2/3)/N, and 63 = q{q — 2{m + m))/N. In Section 
IT^ we take the limit m^m^q N. In this limit, ei = 62 = 1, 63 = 0. 
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Appendix C: Related Models 


2. Other Processes: Kerr, TPA, FCA 


Equations ( 4Qp2 ) are the simplest free-carrier model: 
identical modes, no two-photon absorption, no interac¬ 
tion between carriers, no excitons. In many ways it is 
unrealistic. However, it forms the basis for generalized 
models that include these effects and better approximate 
the real system. 


A host of additional processes may be relevant in 
semiconductor cavities: among the most important are 
the Kerr effect, two-photon absorption (TPA), and free- 
carrier absorption (FCA). Thermal effects and excitonic 
effects, while very important for some systems, are be¬ 
yond the scope of this paper. 


1. Non-Identical Modes 


TPA and Kerr 


The most obvious generalization is to include many 
non-identical carrier modes. This means that, rather 
than grouping all of the modes together into (m,m), 
they are binned into spectrum of modes (ma,ma). A 
similar binning technique is used in many-atom cavity 
QED when the atomic couplings are not equal ^T8\. The 
equations are a straightforward generalization of ( 4Q]|^ ): 


da = 


drris, = 


drris, = 


npT] 


i^c ^ ^ ^ (^a^a T ^a^a) 


adt — d^a 
(Cl) 


|^77a ^ Anr,sJ^Si TrCja^a^a 

+ ^ (Tb^a^b - Ta^b^a) dt + d^^,a (C 2 ) 

b 

l^T^a a a Anr^sJ^aL Arc^aL^aJ^SL 

T ^ ^ ('Tb^a'^b Ta^b'^a) dt + d^m,a 


In indirect-gap materials, like silicon, the linear ab¬ 
sorption is not an effective pathway for carrier genera¬ 
tion. Instead, two photon absorption is the dominant ex¬ 
citation process. Typically, two-photon absorption also 
comes with a dispersive (Kerr) effect. In other cases, 
the band gap is tuned to be very close to the photon en¬ 
ergy, and both processes are important. Unlike linear ab¬ 
sorption, which tends to create carriers very close to the 
band gap, two-photon absorption tends to create highly 
excited carriers with excess kinetic energy. After excita¬ 
tion, these carrier quickly thermalize and subsequently 
decay. 

We can model this with the following Hamiltonian and 
decay process: 


H — ^ Ax(nx + "^x) T W ^ X ^ ^-i-x) 

X 

(C5) 

Ax^k ~ 


The only change here is the introduction of indices 
and the cross-scattering terms 7a^b- These terms, like 
the other carrier excitation / decay terms, have Poisson 
statistics. The Poisson statistics of different modes are, 
of course, correlated just as the flows are - this conserves 
total carrier number in the scattering processes. 

If scattering between bins rua is fast compared to car¬ 
rier excitation or decay, we can replace the mode oc¬ 
cupations by the thermal average rua = fe,ai{T)m,fha = 
fh,ai(T)m^ where /e,a: fh,aL ^rc normalized Boltzmann dis¬ 
tributions. In terms of the total carrier numbers m = 
and fh = we recover Equations (40||42), 

with the effective rates: 


<^ = E'^a/a(T) 

a 

Tnr ~ y^Tnr,a/a(^) 

a 

7rc = Y.^rcMT)UT) 


5 = Y,SMT) 

a 

Tnr ~ ^ ^^ Anr^ajfaiiT^ 
a 

(C4) 


where the new modes Cx, Cx defined for the highly excited 
carriers. Note that, as these modes are highly excited, 
there is no process Tk^x- 

Since the excited state is so short-lived, it can be adia- 
batically eliminated. For on-resonant transitions Ax = 0 
this gives a two-photon absorption term /3; in the off- 
resonant case Ax ^ 0, one finds two photon absorption 
plus a dispersive (Kerr) term. 

These effects add the following terms to the Wigner 
equations: 


A{da) = {—ix — P){oAa)adt — 2 ^/Pa''dPf 3 (C7) 

A(dm) = (i{a*afdt+TP{{oi*f<lPp+(^d(i}){CS) 
A{dfh) = P{a*a)'^dt + (^{a*)‘^d(3 13 + dpp){C9) 


Note how Equations (C7 C9) predict that single 
electron-hole pair is created for every two photons ab¬ 
sorbed, in contrast to linear absorption ( 4Qp2 ), where 
the ratio is one-to-one. 
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b. Free-Carrier Absorption 


In some materials, including silicon, free carriers can 
increase the absorption of the medium, an effect known 
as free-carrier absorption. In addition, for indirect band- 
gap materials, the free-carrier dispersion is larger than 
the band-filling result (46) predicts, due to the collective 
response of the free-carrier plasma [3]. These effects can 
be accounted for by adding the phenomenological terms: 


H = —i(^fcd^ + (CIO) 

— \/Tfca^ ^5 VTfcaT ^ (Cll) 


This can be accommodated in the model (40 42) if the 
substitution Sc + —'<^7fca/2 is made and an extra 
noise is included: 


A{da) = -\/7fcaTO + J(ca.m d(3fci, (C12) 


where d/3fca is another vacuum Wiener process. 

Altogether, the Wigner equations for the free-carrier 
cavity, including two-photon absorption, FCD and 
FCA, take the form: 


da = 


-- - ^(Ac FScmFScfn) 


a dt 


--^adt - ^/ndPr, + {-ix -/3){a*a)adt - 2^/pa*d/3is -^Kca'in + Jica.'m dfi^, 


(^dOi)r] 




dm = [—jnr'mdt — dwm\ + dt + dwrc] “ [a'^{da)rj + a{da)*] 


a'^{da)i3 + a{day^ 


dm = [—+ \—^rc'^'^dt + ^y^c^^ dwrc\ — [<a*(d(a)^ + a{da)'^ — 


a'^{da)p -h a{daY^ 


(C13) 


(C14) 


(C15) 


3. Single-Carrier Approximation 

The single-carrier approximation assumes m = fh = and replaces ScF Sc ^ dc, 7fca + 7fca ^ 7fca- The equations 
reduce to: 


da = 


dNc = 


---i{AcFScNc) 


a dt T 


~adt - y/ndldr^ + {-ix-li){Q*a)adt - 2\/pa*dfip -\/FcF^c dpk, 


(^dcx) r] 


{dcx)^ 


(C16) 

/ - o ! - a"^(da) 3 a(da)% 

[-InrNcdt - YAnrNcdWnr] + [-IrcN; dt + ^dWrc] “ [(a*(da)^ + ^(da)*]- - - - 

(C17) 


Appendix D: Carrier Detuning in terms of Material 
Properties 


der in the optical field the light-matter coupling goes as: 


In this section we derive expressions for the coupling 
constant g\^ and the carrier-dependent detuning dk as a 
function of material properties. This is important be¬ 
cause it allows one to match the results from this work 
to the semiclassical treatment of FCD found elsewhere in 
the literature. 

In standard single-particle electrodynamics, to first or- 


idint = 


eA • p 
mo 


(Dl) 


This can be generalized to many-particle systems 
by “second-quantizing” the Hamiltonian in terms of 
fermionic creation / annihilation operators /k,/k m- 


Hint 


k,lG{states} 


(D2) 
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Consider a two-band model. The /k here represent 
both valence-band and conduction-band states. If the 
field A is driving at optical frequencies, only transitions 
between the valence band and conduction band need be 
considered - for these, the Hamiltonian becomes: 




Too 


+CkCk(k, v\A{x, t) c) (D3) 


For a resonant structure, E{xA) and B{xA) depend on 
the normal-mode fields E^{x) and their time-dependent 
amplitude (t) (which becomes the photon annihilation 
operator when the system is quantized). Working in the 
Coulomb gauge 0(x,t) = 0, £" = —dAjdt and A{xA) is 
given by: 


t) = Re 


-i^/2h/u)eQa^Eu{x)e 


(D4) 


Here E^ is normalized so that J n{x)^\E^\‘^(Px = 1, 
and is the photon annihilation operator. For a good 
resonator, typically only one frequency uj is relevant 
(though multiple frequencies is a simple extension of this 
work), so hereafter we replace a^j ^ a. Going into the in¬ 
teraction picture and neglecting rotating-wave terms and 
adding an arbitrary phase shift to the q to fix the sign 
of ^k, we find: 




= ( 


«'CkCk - acj^cj^) (D5) 


— n A A 


with coupling constant g\^ given by: 


= 


h 


mo V 2ci;eo 


E^{x) • (k,c|^k,i;) 


(D6) 


The electronic and photon parts to the Hamiltonian 
take their canonical forms. The end result is (11). 


Having derived the coupling ^k, we proceed to express 
the carrier-dependent detuning in (46) in terms of ac¬ 


tual material properties. The carrier-dependent detuning 
is what fundamentally limits the performance of a free- 
carrier device - it sets the minimum number of carriers 
needed to switch by one linewidth, the energy figure of 
merit for a photonic switch. It is given by: 


A(m, m) = y^((5kTOk + f^k^k) 


(D7) 


with (5k = 5^/(Ak - Insc), as in (46). 


This section considers two common cases: a HI-V semi¬ 
conductor near the band gap, where band filling is domi¬ 
nant, and silicon far from the band gap, where the plasma 
effect dominates. These effects are well studied in bulk 
materials; the point of this section is to translate them 
to the optical resonator picture used in this paper. 


1. III-V Semiconductor near Band Gap 

Here, the dominant effect comes from band-filling dis¬ 
persion. We assume that all modes have roughly the 
same energy. Eg , and that the optical field is at = xEg , 
where x <1. If x ~ 1, then one can show that the carrier- 
dependent detuning takes the form: 


A(m,m) ^ E 4(mk + mk) 


E 


F e |F/^(xk) ’ Pcv I 


mgeoF’^ 2x(l — x) 


(mk + mk) 


/i^e^|F^(xk)P|Pc 

mleoEl 


\^uj (^k) ' Pcv I 


2x{l — x) 


(mk + mk) 


(D8) 


where mo is the bare electron mass and Pcv is the matrix 
element (k, c| p |k,'z;) between conduction- and valence- 
band states 

This is a two-band calculation, which only includes 
transitions from a single valence band. Adding a sec¬ 
ond valence band doubles the effect of the electrons - 
since each electron “blocks” two transitions, one from 
each valence band, its bandfilling effect is doubled (Fig¬ 
ure]^. This does not happen for holes, since each hole 
only “blocks” the one transition to the conduction band. 
Thus the correct carrier-dependent detuning is: 




mleoEl 


2 mk + mk 


2x{l — x) 


(D9) 


To get a sense of scaling, we replace ^ 

|Fc^P/(^o^)- Here, E^ is designed to have near-unit am¬ 
plitude within the cavity, and V is the mode volume. Un¬ 
like F^, F^ is not normalized (its integral is not one), but 
having near unit-amplitude is what matters here. The 
mode volume is defined in terms of a normalized quan¬ 
tity, V = U/(A/n)^, which is 0(1) for photonic crystals 
and 0(10) for rings. Instead of looking at A, we look 
at A/cj, since this is unitless, and we are well aware that 
A(m, fh)/uj ^ 1/Q means that enough carriers have been 
injected to move the cavity one linewidth. 


A(m, m) 

UJ 


3 e‘^no\pcv\‘^ X 

Sn^m^hc^eoV 2(1 — x) 


E 


|F(xk) • pcv I 


2 mk + mk 

3 


(DIO) 

This is a product of four terms. (1) The first is a 
bunch of fundamental constants, plus material and cav¬ 
ity parameters like the cavity size and the index of re¬ 
fraction, and the magnitude of the matrix element \pcv \ • 
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FIG. 15: Bandfilling in a direct-gap III-V semiconductor. 
Carriers block certain optical transitions, changing the ab¬ 
sorption spectrum, which in turn alters the index of refrac¬ 
tion. 


derivations of the band-filling dispersion EES], under the 
replacements \pcv\‘^ Egm^l2mQ and {\E{x)\‘^)\^ 1 

(this is always 0(1) and the equality can be imposed by 
scaling V). Because of the rotating-wave approximation 
taken in this paper, it will not be valid when x deviates 
far from 1. However, optimized devices exploiting band¬ 
filling always operate near the band gap. 


2. FCD in Silicon 

In silicon, the indirect band gap makes the band-filling 
effect very weak. Instead, free-carrier dispersion is dom¬ 
inated by the plasma effect [3] . Consider a simple Drude 
model with a carrier density N. The index of refraction 
is modified as follows: 


These constants depend only on the device, not on the 
wavelength operated at or the particular carrier mode 
being excited. (2) Next we have a term that depends 
on the closeness to the band edge: x/2{l — x). In prac¬ 
tice, we will want x to be as close to one as possible in 
order to maximize this quantity - however, as x ^ 1 lin¬ 
ear absorption limits the cavity Q, so there is a tradeoff 
here. (3) Next is a position term that depends on the 
field strength at Xk, the location of the carrier (assuming 
carriers localized to well below a wavelength here). (4) 
Finally, the carrier number. 

When carrier thermalization and diffusion is fast com¬ 
pared to the decay processes, one can write this as an 
effective single-mode model, where the per-carrier detun¬ 
ing is given by the thermal average: 


2 , 2 
n ^ Uq 



nce^/wcn^eo 
iujIT 


(D12) 


where ric = n^p are the densities and rUc = rUg, rrih are 
the masses for electrons and holes. In the high-frequency 
limit uj ^ Up, the real part dominates and this becomes: 


2noeoE‘^ 


n 

me 


P 

rrih 


(D13) 


Assuming the carriers are confined to a volume V, and 
defining the dimensionless V = V/{X/n)^ as above, and 
using standard coupled-mode theory to convert An to a 
detuning, we find: 


A(m, m) 


u 


Se^nolpc 


S'K'^m'khc^eoV 2(1 — x) 


\E{x]i) -pcvl 


2m + m 


e^nolpc 


X 2m-\-m 

^ STT^mlhc^eoVW^ 3 ^ ^ 

In the limit x ^ 1, this is consistent with previous 


A(m,m) = 


e^no 


Wn^hc^eoV 


m m 

- 1 - 

rUe mh 


(D14) 


A more detailed treatment shows that the dependence 
is linear for electrons, but nonlinear for holes [34]. This 


nonlinearity can be treated phenomenologically in (40 


42); the quantum noise terms derived in this section do 
not change. 
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